LINEAR VS. NONLINEAR EFFECTS FOR NONLINEAR 
SCHRODINGER EQUATIONS WITH POTENTIAL 



REMI CARLES 

Abstract. We review some recent results on nonlinear Schrodinger equations 
with potential, with emphasis on the case where the potential is a second order 
polynomial, for which the interaction between the linear dynamics caused by 
the potential, and the nonlinear effects, can be described quite precisely. This 
includes semi-classical regimes, as well as finite time blow-up and scattering is- 
sues. We present the tools used for these problems, as well as their limitations, 
and outline the arguments of the proofs. 



1. Introduction 

This paper is a survey of some recent results on nonlinear Schrodinger equations 
with potential. A particular attention is paid to the case where the potential is a 
second order polynomial. In this case, the fundamental solution is known explic- 
itly, through a generalized Mehler's formula; the linear dynamics is well understood. 
The most important remark is that we can assess the action of some Heisenberg 
observables (which can be exactly computed) on a class of nonlinearities. With this 
tool, we can understand the interaction between the linear effects caused by the po- 
tential, and the nonlinear effects. In the semi-classical regime e — > 0, we emphasize 
some critical scales measuring this interaction, and describe the critical phenomena. 
When e = 1, we observe the effect of such potentials on finite time blow-up; the 
potential may "create" some blow-up (think of the harmonic potential), or delay, 
even prevent, this phenomenon ("repulsive" harmonic potential). 

We give two motivations to study nonlinear Schrodinger equations with potential. 
The first one arises from physics, for Bosc-Einstein condensates. The model involves 
the equation 

(1.1) ihd t u h + ^h 2 Au h = V{x)u h + ah 2 \u n \ 2a u h ; x £ R n , a e R , a e N , 

where the role of the potential V is to confine particles. The cases most currently 
considered are when V is quadratic (isotropic or anisotropic harmonic potential), 
when V is lattice periodic, or when V is the sum of two such potentials (see e.g. [3H1 
ElEZIEj)- We shall not discuss the physical relevance of this model, but notice that 
explicit mathematical formulae are available in the case of the harmonic potential. 
Our main motivation is more from a mathematical point of view. The presence of a 
potential alters the propagation of the wave in the linear case. Similarly, nonlinear 
problems may lead to typical phenomena, such as finite time blow-up. How can 
these two effects interact? This vague question seems to raise many complex issues. 
We give some very partial answers, essentially restricted to the case where the 
potential is a second order polynomial. These examples may be viewed as a first 
step for a general study supporting or contradicting the intuition. This may be 
compared with the approach of N. Burq, P. Gerard and N. Tzvetkov, who analyze 
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the role of geometry for nonlinear Schrodinger equations on compact manifolds (see 

e.g. HE1). 

The initial value problem we study is 

(1.2) ied t u s + ie 2 Au £ = V{x)u e + A £ |u £ | 2<T u e ; u 8 |, n = ug, 

2 it— u 

where e G]0, 1], x £ R n , V" G C°° (M™ ; K) , A e G R, and cr > with a < ^ if n > 3 
(the nonlinearity is iJ 1 sub-critical). We consider two regimes for the parameter s: 

• The semi-classical limit e — * 0. This regime gives hints to understand high- 
frequency phenomena and provides us with tools which are extremely useful 
in the nonlinear case, when V is a second order polynomial. 

• The case e = 1. This case is better understood thanks to the semi-classical 
analysis. We study in particular the global existence issue and two of its 
companions: finite time blow-up and scattering. 

We shall not discuss any regularity issue here, and always assume Mq € S, where 

S := {/ G S'(R n ) ; ||/|| s := \\f\\ L 2 {Rn) + ||V/|| w) + ||x/|| i2(K n) < +00} . 
We do not discuss the question of solitons either. We denote 

Hy = + V(x) ; U v {t) = e^ H v , and simply H v and U v when e = 1. 

This paper is organized as follows. Section |21 is devoted to general results. First, we 
recall some classical results in the case V = and e = 1 , as well as some techniques 
to prove them. We then discuss which ones can be easily generalized when V is 
not identically zero, and present some cases where a change of variables makes it 
possible to relate the case V = to the case where V is not trivial. 

In Section |2| we recall the generalized Mehler's formula, and motivate the in- 
troduction of some particular Heisenberg observables, from a nonlinear point of 
view. 

In Sectional we state some results about (|1.2|l in the semi-classical limit; we 
outline the techniques, and discuss their limitations. 

In Sections El and El we set e = 1, and analyze the influence of the potential V 
on finite time blow-up and global existence issues. 

2. General setting and consequences 

2.1. Some results on the nonlinear Schrodinger equation. In this paragraph, 
we assume e = 1, and recall a few things about the Cauchy problem 

(2.1) id t u+^Au = X\u\ 2a u ; u\ t=0 =u , 

where x G R n , A G R, and a > with a < ^2 if ^ > 3. All the results we mention 
can be found in |14) . 

Theorem 2.1. Suppose uq G H 1 (M n ), and a > as above. Then there exist 
T*,T* > 0, and a unique maximal solution 

u G C (] - T*, T*[; H 1 ) n L« oc (] - T„ T*[; W^ 2a+2 ) 

to 112. lj l. where q — 4 ^ 4 . It is maximal in the sense that ifT* < 00, then 

lira \\V x u(t)\\ L 2 = +00 . 

In addition, the following quantities are independent of time: 
Mass: \\u(t)\\ L2 = \\u \\ L2 , 

^2.2^ ]_ 2 /\ 2cr I 2 

Energy: E a := - ||V K u(t)|| L2 + — — ||it(t)|| L2o -+2 = amsf. 
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This result was first proved in [23] (with a slightly different statement), and 
revisited in |26l 155) . The proof relies on a fixed point argument on Duhamcl's 
formula 

(2.3) u(t) = U (t)u - iX [ U {t~s)(\u\ 2a u){s)ds, 

Jo 

where following the notation introduced in Section Uo(t) = e l i A . The modern 
tool to prove this result is Strichartz estimates, after j^- We recall a statement for 
the case of Schrodinger equations (551 134| (see ^U] for the following adaptation): 

Lemma 2.2. Let (?7(t)) tgK be a unitary group on L 2 (R"), satisfying the dispersive 
estimate 

(2.4) \\U(t)\\ L i_ LOO < w(i)™ /2 , with w > and u e Ll{R) . 
We say that a pair (q, r) is admissible if 

2 n n 2n 
— I = — : q, r > 2, with r < 



q r 2 n — 2 

Then for any T G R + , and any admissible pairs (q,r) and (q,r), we have 

\\U(t)f\\ Lq{] _ TT[ . Lr) < ||wl]_2T,2T[||^ ||/|| L 2 , 



t 

U{t-s)F{s)ds 



1/9+1/5 

L 1 N*' WL'i' (]~T,T[;Lr') ■ 

Li(]-T,T[;L r ) 



< wl 1 _ 2T , 2r r||;?^ /, '||F|| 



It is straightforward that Uo(t) satisfies the assumptions of Lemma 12.21 with 
w(t) = (27r |i | ) 1 , from the formula 

(2.5) Uomx) = -^J^^ mdy . 

Theorem 12.11 then follows from Strichartz and Holder's inequalities, as well as 
Gagliardo-Nircnbcrg inequalities when a > ^. A similar result holds when «o € S; 
essentially, we can remember that in addition, u € C(R; E). 

In the H 1 or E case, global existence can be deduced from the conservations 
1)2.2(1 . when A > (defocusing nonlinearity) , or A < and a < ^ for instance. If 
A < and <j > — , finite time blow-up may occur: 

Proposition 2.3 (Virial Theorem, |27]). Assume A < and a > ^ (with a < — ^ 
if n> 3). If uq S E is smc/i i/iai i?o < 0, i/ien solutions to <|2.1|) 6Zow up in finite 
time, in the future and in the past: T* and T* are finite. 

The idea of the proof for this result is to introduce the function y(t) = \\xu(t)\\^ 2 
and to notice that under the above assumptions, y(t) < Eq < 0. Since y(t) > 0, 
this proves that a singularity appears, both for t > and for £ < 0. 

A slightly different approach to recover this result is to use the Galilean operator 
J(t) — x + i£V X ) an d the pseudo-conformal conservation law, discovered in |25| : 

(2 ' 6) Jt Q 11 J( * H " 2 + ^TT* 2 |I U (*)H^) = ^Tl( 2 - na)t\\u{t)\\%i%. 

As noticed in [51], expanding the above formula, one retrieves Proposition 12.31 
Another application of the pseudo-conformal conservation law, as motivated in 
J25I, concerns scattering theory: 

Proposition 2.4 (Scattering theory in E, [2312311211301 ) - Assume < a < 
if n> 2, a > 1 if n = 1. Assume in addition A > 0. Denote 



2-n + y/n 2 + 12n + 4 
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• Let u_ € E. There exists a unique u G C(R; E) solution of (|2.1|l . smc/i i/iai 

lim || t/" (— - "-lis = 0. 

i — > — oo 

• Let tio e E. If a > oo(n) or £/ ||uo||e *s sufficiently small, then there exists 
a unique u+ G E swc/i that if u € C(K; E) is t/ie solution of ((2.1(1 . i/ien 

lim ||tT (-t)«(t) - u+|| s = 0. 

i — >+oo 

We denote by S : u_ i— > w + £/ie scattering operator. 

We emphasize two properties of the Galilean operators, from which the above 
results (Equation 1(2.6(1 and Proposition ^. 4(1 follow: 

(i) J(t) is an Heisenberg observable: 

(2.7) J(t) = U {t)xU (-t). 

(ii) It also reads 

(2.8) J(t) =ite ii &V x (e- iS # 
As a consequence of 1(2.7(1 . J(t) commutes with the linear part of l|2.1|) . 



(2.9) 



J(t),id t + ^A 



0. 



The second point has two interesting straightforward consequences: 

(ii)' Weighted Gagliardo-Nirenberg inequalities: for 2 < p < (2 < p < oo 
if n = 1), there exists C p depending only on n and p such that 

(2-10) ll/ll w < ll/ll^) || JW/Hg^) ; *(p) := n Q - ^) • 

(ii)" If F(z) = G(\z\ 2 )z is C 1 , then J(t) acts like a derivative on F(w): 

(2.11) J(i) = d z F{w)J(t)w - i\F{w)J{t)w . 

Roughly speaking, ((2.9(1 and ((2.11(1 make it possible to have estimates for J(t)u 
in the same way as for V x u, when u solves (12.1(1 . Then ((2.10(1 yields dispersive 
estimates for the nonlinear equation, which are the same as in the linear case, 
provided that J(t)u G L^°L 2 . These arguments are the key ingredients for the first 
point of Proposition l2~4l (existence of wave operators). Then ((2.6(1 yields estimates 
on || J(i)u||x,2 which, along with the conservation of mass ((2.2(1 . prove the second 
point of Proposition ^. 41 (asymptotic completeness). 

Remark. The properties ((2.911 . 12.10(1 and ((2.11(1 are the analog of those satisfied by 
the conformal Killing vector-fields used to study the wave equation (|36|'l. 

2.2. Introducing a potential. We now turn to the case of 1(1.2(1 . Drawing a 
parallel with the case of 1(2.1(1 . the first tool we seek is Strichartz estimates for Uy- 
This has been, and this is still, a very active area of research. The possible presence 
of eigenvalues shows that in general, one cannot expect not only the same global 
dispersion as for [/o, but also global in time Strichartz estimates. In a more subtle 
way, resonances may also be an obstacle to dispersion ([SOI)- Many results have 
been obtained though, and we refer to the introduction of |48( for a very nice review. 
In (48j . Strichartz estimates are obtained for potentials which may depend of time, 
a case we do not consider. Notice that these involved results rely on perturbation 
arguments (see @2| f° r an interesting exception). We would like to consider the 
case where the potential may really change the dynamics of the Laplacian, and as 
introduced in we assume that V G C°°(R™;M). 
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It is well known that if V(x) > — a|x| 2 — b for some a, 6 > 0, then Hy is essentially 
self-adjoint on C^°(R™), and that this is a sharp result: if V(x) — ~\x\ 4 for instance, 
this property fails (classical trajectories can reach infinite speed, see ^ElEOJ)- On 
the other hand, if V is positive and grow faster than quadratically, then at least 
for n — 1, the kernel of Uy is nowhere C 1 f |56p. but smoothing properties make it 
possible to solve the nonlinear problem (|1.2() in some cases f |57|). We shall restrict 
our attention to sub-quadratic potentials, for which it is possible to construct a 
parametrix, locally in time: 

Proposition 2.5 ((22E3)- Let V £ C°°(R n ;R) be such that d a V £ L°°(R") for 
\a\ > 2. Then there exists 8 > independent of e s]0, 1] such that for \t\ < 5, 

(2.12) U v (t)f(x) = e- in ? B * at \ n/2 [ k £ (t,x,y)e^ t ^/ £ f(y)dy, 

\2TTet\ 1 JR" 

where S solves the eikonal equation 

(2.13) d t S+^\V x S\ 2 + V(x) = 0, 

and d%d^k £ £ L°°(] - 6,S[xR 2n ) for all a, (3 € N n , uniformly in e G]0,1]. 

This result relies on perturbation arguments: for small times, the influence of 
the potential on classical trajectories is controlled, hence a formula similar to H2.5|l . 

Since we know that Uy is a unitary group on L 2 , the above result shows that 
for \t\ < 5, U v is dispersive: 

IIW)llL^L-^pp^- 

We infer that there are local in time Strichartz estimates. Notice that in general, 
one cannot expect global in time estimates: when V is the harmonic potential or 
is lattice periodic, Hy (hence H v ) has eigenvalues. 

One can then mimic the proof of Theorem 12. II Unlike for the case of ^2.1|l . the 
gradient does not commute with the linear part of the equation: 



ied t + ie 2 A - V(x) ) V x u £ = u e W x V + X E W X (| 



u e \ 2,T u £ ) 



The new term is u e \7 x V. Since V is sub-quadratic, |VV(x)| < 1 + |a;|, which 
suggests to consider xu £ as a third unknown, after u £ and \7 x u e . We then have 

(ied t + \e 2 A - V{x)^j xu e = V x u e + \ e x (\u E \ 2a u E ) . 

It is easy to prove: 

Lemma 2.6. Suppose that V £ C°°(R™;R) is sub- quadratic, Uq £ E, X £ £ R, and 
a > with a < — ^ ifn>3. Then there exist T £ > (depending on n, a, V, \X £ \ 
|| u ollsJ; and a unique solution 



u £ £(7(]-r,f[;S)nL? oc (] — T £ , T £ [. W 1 ' 2a+2 ) 
to l|1.2|l . where q — 4 ^ 4 . The following quantities are independent of time: 
Mass: \\u £ (t)\\ L2 = \\u £ \\ L2 , 

Energy: E £ v := \ \\eV x u £ (t)\\ 2 L2 + \\u £ (t)\\ 2 L a ++ 2 + [ V(x)\u £ (t, x)\ 2 dx . 

Z (7+1 J Rn 

In addition, if a < ^ (L 2 sub-critical case), then one can take T £ = +oo. 
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Notice that without further assumption on V, we cannot state a criterion for 
the obstruction to global existence, as in Theorem 12. II Global existence in the L 2 
sub-critical case follows from the same ideas as in [53]: one has global existence at 
the L 2 level, from Strichartz inequalities (local in time inequalities are sufficient) 
and the conservation of mass. One deduces global existence in E by considering 
|u £ | 2cr as a potential in the equations for V x u E and xu e . 

If in addition V is non-negative, then the assumption Uq S S can be weakened: 
it suffices to consider initial data in the domain of \J — A + V (see In that 

case, the obstruction to global existence is the same as in Theorem 12. II 

2.3. Removing the potential. It turns out that for some specific potentials, an 
explicit change of variables makes it possible to relate the solutions of (|1.2I) to the 
solutions of the same equation with V = 0. These potentials are: 

• The linear potential (Stark potential): V(x) — E ■ x. 

• The isotropic harmonic potential V(x) = \x\ 2 , in the case where the non- 
linearity is L 2 -critical, a = — . 

• The isotropic repulsive harmonic potential V{x) — — \x\ 2 , in the case where 
the nonlinearity is L 2 -critical. 

Introduce v s the solution to the equation with no potential: 

ied t v s + ±E 2 Av e = \ s \v e \ 2 °v e ; « £ | t=0 ="o (=« £ U) ■ 

2.3.1. Linear potential. Assume that V(x) = E-x, where E £ R™ is constant. Then 
as noticed in ^3], the Avron-Herbst formula, discovered in the linear case X s — 
(PP), relates u e and v e : 

(2.14) u%t,x) = v s (t,x + ^E^ e -*(«*-+£l*l>. 

Therefore, one can use Theorem 12 .11 to deduce local existence results in H 1 and see 
that the only obstruction to global existence is the unboundedness of V x u £ in L 2 . 
Similarly, the Stark potential does not change the possible blow-up time, nor the 
scattering theory (replace Uq by Uy in Proposition 12.4(1 . In the case of finite time 
blow-up, the linear potential only shifts the set where this phenomenon occurs: if 
v s blows up on a set X at time T e > (see e.g. [3B] for a definition of such a set), 
then u e blows up (at time T e ) on X — ^ T J E. 

2.3.2. Isotropic harmonic potential and conformal nonlinearity. Assume that a = ^ 
and V(x) = \ \x\ 2 with ui > 0. Then provided that the right hand side is defined, 

1 _-,#L^ tan(J ( , fttmujt x 



(2.15) u £ (t,x) = - — e~ l ^ x tanw V . 

y J y ' (coswi)™/ 2 V w cosu)t / 

This was first noticed in for the linear case (A £ = 0), and in El fo r the 
nonlinear case with critical nonlinearity. Note that despite this expression, u £ 
needs not blow-up at time t% = (when the cosine is zero). Indeed, as t — > t\, the 
time variable tanu,t for v £ goes to infinity. Then dispersive properties of solutions 
to 1)2.1(1 may compensate the cancellation of the cosine. 

However, it is true that the harmonic potential generates more blowing-up solu- 
tions (the criterion for finite time blow-up is the same as in Theorem 12 .1|) . First, if 
v e blows up at time T e > 0, then u e blows up at time arcta ^(" T 1 ; w hich is always 
smaller than T 6 . Second, if v e is not a dispersive solution, typically of solitary wave 
v e (t,x) = e JKt / £ Q(-), then u e blows up at time j-. From a heuristic point of view, 
the confining properties of the harmonic potential are sufficient to concentrate an 
energy which is not naturally dispersed. 
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Similarly, because the harmonic potential prevents the solutions from being dis- 
persive as in the case V = 0, no scattering theory must be expected. 

2.3.3. Isotropic repulsive harmonic potential and conformal nonlinearity. Assume 
that a = - and V{x) = — ^\x\ 2 with u > 0. Then as noticed in [S], 

(2 16) uHt x) = - e ^l*l 2 ta„hM V ( tanhM)^ x 

(cosh(^)) n/2 V w cosh(wi) 

A criterion for global existence is not obvious in this case, since the potential is 
unbounded from below. It was proved in |H] that for an isotropic repulsive harmonic 
potential (and a general nonlinearity as in l|1.2(l ). the obstruction to global existence 
is the same as in Theorem 12.11 (see also Section EJl . Opposite to the harmonic 
potential mentioned above, the repulsive harmonic potential tends to prevent blow- 
up. If v £ blows up at time T e > 0, then u £ blows up at time algta °^" T ) if w T e < 1. 
That means that for "small" values of uj, blow-up is delayed. On the other hand, 
if uiT £ > 1 , then u e does not "see" the blow up of v e , and is global. If wT e > 1 , u e 
has even exponential decay as t — > +oo. The limiting case wT E = 1 is interesting: 
it is known that there exists no (L 2 ) solitary wave ([SU E3 EH])- The analog of a 
solitary wave is a solution that grows exponentially in time (see [HJ). 

Unlike the harmonic potential, the repulsive harmonic potential tends to spread 
out a solution to 12.1f> : we saw how it changes the finite time blow-up phenomenon. 
It is not surprising that things go very well as far as scattering is concerned, thanks 
to exponential decay. We study this issue more precisely in Sectional 



3. Mehler's formula and applications 

When V is a second order polynomial, a lot can be said. First, the linear case 
is extremely favorable, because U v is known explicitly: this is the (generalized) 
Mehler's formula f |21M29| 1. If V(x) = ^ OjkXjXk + ^i x j + c j then reducing the 
quadratic part, a change of orthonormal basis (which leaves the Laplace operator 
invariant) and a change of the origin simplify the expression of V: 

V{x) = ^2 ( 5 i~2 X 'j + h J x 3 ) + c ; Sj € {-1, 0, 1}, ujj > 0, Sjbj = . 

Using Avron-Herbst formula Ij2.14|l . we can get rid of linear terms bjXj. Taking 
e ict/e u e ag a new un k nown function also removes the constant term. We therefore 
assume that V is of the form 

" uj 2 

(3.1) V(x) = £ S^x 2 ; 6j e {-1, 0, 1}, w s >0. 

i=i 

The Mehler's formula then reads: 

n , v 1/2 r 

(3.2) U v {t)f{x) = JJ / e^y^f(y)dy, 
where 



™ 1 / 

s{t ' x>y) = ^W){ 



X 2 + y 2 \ 
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and the functions gj and hj, related to the classical trajectories, are given by: 

' smh(ujjt) 



(3.3) (g 3 (t),h 3 (t)) = { 



sin(wji) 



UJjtU 


, iffy 


= -1, 


(M) 


, iffy 


= 0, 




, iffy 


= +1. 



Recall that if there exists fy = +1, then e~ ltH v has some singularities, periodically 
in time (see e.g. [H]). This affects the above formula with phase factors we did 
not write (which can be incorporated in the definition of (^ (i)) 1 / 2 ). 

Let us examine the Strichartz estimates we can deduce. If fy = +1 for any j, then 
only local in time Strichartz estimates are available. This is not surprising, since the 
harmonic oscillator has eigenvalues. On the other hand, if ujj = —1 for any j, then 
we have ||E7y-(t)||£i_ > £oo < |27r£i| - ™/ 2 , an estimate which is independent of all the 
Uj's: we have the same dispersion as in the case with no potential (recall from the 
previous section that the repulsive harmonic potential accelerates the "particle"), 
and global in time Strichartz estimates follow. Actually, this property remains if 
fy = —1 for at least one j. This stems from the formulation of Lemma 12.21 the 
function v(t) given by Mehler's formula is in the weak L 1 space for small values of 
\t\, and in the strong L 1 space for "large" values of \t\. 

Remark. We would like to point out that in the case where fy = — 1 for at least one 
j, we have global in time Strichartz estimates, while there are trapped trajectories. 
The one-dimensional case shows the mechanism. If n = 1 and to = 1, the classical 
symbol of H v is ^(£ 2 — x 2 ). Computing the Hamilton flow, we have 

(+\ u 4. , c ■ v, + x ° + £° t , - £o - t 
x(t) = xq cosht + £0 sinht = — - — e H — e 

If 20 + £0 = 0, then the trajectory is trapped in the future. This aspect can be 
compared to the results of ^3 > where it is proved that smoothing effects occur 
in the future provided that the classical trajectories are not trapped in the past. 
However, the results of |16| include potentials which grow at most linearly in x, 
and JB] does not consider the case of potentials. On the other hand, smoothing 
effects yield another approach to prove Strichartz estimates (see e.g. jHOj)- In 
our case, there are trajectories trapped in the past or in the future, but global in 
time Strichartz estimates are available. It seems that the link between classical 
trajectories and (global in time) Strichartz estimates remains to be clarified. 

To introduce the tools we use in the rest of the analysis, recall the approach 
followed in [S]. We consider the case of an isotropic harmonic potential, V(x) = 
^rr|a;| 2 . In the linear case X £ — 0, u £ is therefore given by Mehler's formula. A 
formal stationary argument (which can be justified) shows that if Uq = f does not 
depend on e, then for \t\ < ^j, 

In the nonlinear case, we expect that if, say X s = 0(e fc ) with k sufficiently large, 
then nonlinear effects should be negligible in the semi-classical limit, at least before 
the first singularity at time This is proved in [S], and we recall the argument 
in Section ^ A natural candidate for an approximate solution is then given by 
To prove the approximation is valid in the nonlinear case, L 2 estimates are 



u 



app- 



not sufficient: L p estimates are needed, for other values of p, and one can think of 
Gagliardo-Nirenberg inequalities, \\g\\Lp ^ IMIis 11^011*2 • 
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However, these inequalities yield ||w| pp (i)||Lp ^ e~ s ^ p \ which is terrible in the 
limit £ — > for \t\ < This "bad" power of e stems from the highly oscillatory 
phase. Note that the nonlinearity we consider, F(z) = \z\ 2<T z, does not create new 
harmonics in a single phase WKB asymptotics, if only one harmonic is present 
initially. Therefore, in our case, there is only one phase and one harmonic to 
take care about. This suggests to replace the gradient in Gagliardo-Nirenberg 

inequalities by the operator ( e luJ ^~ tan ( wt ) .J. In view of the expression for 
w| pp , we introduce more precisely: 

J £ (t) = icos(^)e- Ml ^ tan ^V x ^^ tan(wt) = -u 2 sin(ujt)~+ico$(ujt)V x . 
Introduce the operator 

r . K , sin(u)t) 
H £ (t) = cos(ujt)x + i — - — -eV x . 

We can then rewrite the energy Ey (which is constant from Lemma EIJ), &s 

E £ v = \ ||eJWH 2 L2 + %- Wmfv + IK(*)ll£& ■ 

More generally, if the potential if of the form 13.1(1 . we define: 

We set J £ = (Jj)i<j< n and H £ = (-ff|)i<j<n- These operators have been known 
for years, in the linear theory (see e.g. (521 147j ). They are Heisenberg observables: 

J £ {t) = U £ v {t)iV x U £ v {-t) ; H £ {t) = U £ ,{t)xU £ v {-t). 

The fact that such Heisenberg observables can be computed exactly is due to the 
assumption that the potential is a second order polynomial. Note that if Sj = 0, 
then we recover the two operators introduced in Section l2~Tl the gradient and the 
Galilean operator. We will see in the next section that when V is of the form 1(3. ljl . 
then J £ and H £ satisfy the same properties as those emphasized in Section l2~Tl 

4. Semi-classical analysis in a nonlinear framework 

The main features of the operators J e and H £ are the following: 

Lemma 4.1. Let V be of the form 1(3. The operators J e and H £ satisfy: 
(i) They commute with the linear part of ((1.2|) : 



J £ (t),ied t + ^s 2 A~V(x) 



H £ {t),ied t + ^e 2 A-V(x) 



0. 



(ii) They can be written as 

H £ (t) = ig j (t)e** a <- t ' x V e d j ( e -**»(*.*)/= .) , 

for some real-valued phases <f>\ and 4>2, solutions of the eikonal equation 1(2.13(1 . 
(ii)' They yield weighted Gagliardo-Nirenberg inequalities, such as 

n 

ii/iu* ^(nifeiWi^ii/iiii^ikw/ii^. 

3=1 

(ii)" IfF(z) = G(\z\ 2 )z is C 1 , then J(t)F(w) and H(t)F(w) are given by lj2~TTJ) . 



10 



Ft. CARLES 



The first point is trivial, since J e and H e are Heisenberg observables. The last 
two points follow from the second one, which is the way we found the operator J e 
in the case of the isotropic harmonic potential. We can consider that two algebraic 
miracles occur: first, we can compute explicitly some interesting Heisenberg ob- 
servables. Second, these operators can be written like (ii), which is a nice property 
in view of nonlinear problems. We discuss these aspects further into details below. 

We now come to the issue of critical scales in the semi-classical analysis of l|1.2Jl . 
We describe two cases: 

• V is an isotropic harmonic potential, and Uq = f does not depend on e. 

• V is of the form (|3.1|l and Uq is a concentrating profile. 

The first case corresponds to the one that led us to introduce the operators J s and 
H e . In the linear case, due to the harmonic potential, the solution focuses at the 
origin at time t = . More precisely, the geometry of the propagation in the limit 
e — > is given by the Hamilton flow. The classical Hamiltonian in this case is 
p(t,x,r,£) — ^{\£, 2 \ + w 2 |x| 2 ). The classical trajectories (rays of geometric optics) 
are given by 

, \ i \ j- sinfwt) 
x(t) = xo cos(wt) + £o • 

LU 

Since / does not depend on e, there is no initial oscillation: £o = 0, and rays meet 
at the origin, periodically in time. When X £ = 0, one has the sharp estimate: 

S(p) 



1 



I cos(wi) I + e | sin(wi) | _ 

This follows easily from Lemma f4. II and the conservation of the L 2 for solution to 
linear Schrodinger equations. Now assume that X 6 — e a , for some a > 0. Then the 
conservations of mass and energy ( Lemma 12. show that the solution u e is global 
in time. For a large, the linear solution is expected to be a good approximation 
for the nonlinear solution. More precisely, two regimes must be considered. Before 
focusing, the solution is of order 0(1), and one can apply WKB methods. This 
leads to a linear approximation if a > 1, nonlinear otherwise. Next, at the focus, 
the linear solution is of order e~™/ 2 , more precisely, we infer from l|3.2[l : 

U ^h X ) eZAlY' 2 Kt) ' where /(0=^ n e--«/W^. 

The nonlinear term can be viewed as a potential £ Q |u e | 2<T . Plugging the above 
asymptotics suggests that a critical value for a at the focus is a = no. 

With the arguments given in it is possible to prove that if a > max(l, ncr), 
then the nonlinearity is negligible in the limit s — > 0, locally uniformly in time. The 
case a = no > 1 is shown to be critical: 

Theorem 4.2 ( 9 ). Let V(x) = and 2 < r < with r = oo if n = 1. 

Assume that the nonlinearity z i— > |z| 2<T z is twice differentiable, and that X £ — e na , 
with no > 1. Assume moreover that o > oo(n). Then for k G N, the asymptotics 
of u £ for \ + (k — l)ir < t < -| + kir is given, in L 2 n L r , by: 



+0 (27r| cost|)"/ 2 V cost / 

where S k denotes the k-th iterate of S (which is well defined under our assumptions 
on f), and <t/j_(x) := {2iir)- n / 2 f{x). 

Note that the phase shift — n^, appearing at each focus crossing, is a linear 
phenomenon (Maslov-Keller index, see ^1]). Thus the only nonlinear effects at 
leading order occur at the focuses, and are described, in average, by the scattering 
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operator S associated to (|2.1|l . We give the main ideas of the proof. The first 
step consists in proving that before the first focus at time t = S, the function 
m| pp defined by (|3.4|) is a good approximation for the nonlinear solution u £ . This 
is equivalent to justifying a WKB asymptotics at leading order, and relies on the 
sharp estimate given by the operator J £ and the weighted Gagliardo-Nirenberg 
inequality it provides (point (ii)' in Lemma |4.1|1 . We can prove that u| pp remains 
close to u £ (in a space which is essentially E), up to time ? — Ae, in the limit 
A — * +00, that is, until focusing effects become relevant at leading order. This 
shows that as predicted, the assumption a > 1 makes the nonlinearity negligible 
outside the focus. 

When t — \ = 0(e), u £ tends to concentrate near the origin, at scale e; V{x)u £ 
becomes negligible, while the nonlinear potential e na \u £ \ 2 ' T is of order 0(1). This 
can be proved thanks essentially to the operator H £ . Note that for t — ? = 0(e), 
the operator J £ and H £ can be replaced respectively by | +i (t — -|) V x and eV x , 
which is another way to check that the potential is negligible. The scaling 

(«) ■ F < f -">-i>(4 i 7) 

turns the description of the caustic crossing into a continuity issue for a scattering 
problem. We have ||£/o(— *)(V' e (0 — ^ , (*))l|i 00 (R;S) = 0(1)7 where tjj solves Ij2.1|l . with 
as a Cauchy datum at t = —00. This scattering state stems from the transition 
between the two regimes for u £ , which occurs for t = ^ — Ae, in the limit A — * +00. 

Next, u £ leaves the focus in a way described by the second point of Proposition ^. 41 
and l]4.1[l . The analysis is symmetric to the one before the focus; we prove the 
asymptotics until t = 7r, where the situation is similar to that at time t = 0. 
Iterating the analysis yields Theorem l4.2l 

Remark. With no additional effort, one can treat the case of initial plane oscilla- 
tions. Let £0 G R™, and introduce 

u £ {t,x) = u £ (t,x - ZoSmt)e< x -^ sint )-Z° cost/£ . 
1 1 2 

If u £ solves (|1.2f) with V(x) — ^j±- and Uq — f, then u e solves (ll.2|l with initial data 

The second case we analyze is when V is of the general form (|3.1() and 
(4.2) u £ (x) = -±p<p (J) (<p € E) ; A £ =e" CT . 

If A e = e Q with a > max(l, ner), then one can show that the evolution of u £ is 
linear, at leading order. The case A £ = — £ n<J has been studied by several authors 
(see e.g. EE!); m tne case a < f an d = — a;o) eia: ^ , where R is a 
ground state. In that situation, u £ evolves as a concentrating profile, with profile 
R, along the Hamilton flow associated to Hy with data (xo,£o)- Heuristically, this 
is so because there is a balance between the dispersive effects associated to Hy, 
and the nonlinear effects (in the case X £ < 0, the nonlinearity is attractive). 

In the case A e = +e™ cr , the two effects mentioned above tend to cumulate, and 
dispersion alters the shape of u £ . Note that even though A e > 0, global existence 
for u £ is not obvious, since V is not necessarily signed. 

Theorem 4.3 Let V of the form (TTT|l . and Uq given by Q4.2J1 ■ Assume 

either that there exists j such that Sj ^ 1 , or that Sj — 1 for all j and the loj 's 
are not pairwise rationally dependent. Suppose that the nonlinearity z 1— > \z\ 2a z is 
twice differentiate. Then the following holds. 

1. For any T > 0, there exists e(T) > such that for < e < e(T), (|1.2|) has a 
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unique solution u s £ C([—T,T]; £). 

2. This solution satisfies the following asymptotics. 

• For any A > 0, 



\L 2 



limsup sup ( \\u E (t) - v e (t)\\ L2 + \\eV x u E (t) - eV^ir (t)^ 

e^O |i|<Ae V 

+ ||f« e (*)-f« e (*)|| J . a )=0 ) where v*(t,x) = 0, X - ) 

and t/> € C(R; E) is i/ie solution to l|2.1|l suc/i £/iai ^|t = o = 
• Beyond this boundary layer, we have 

limsup sup ( ||n £ (i) — it±(*)|| r2 + \\£V x u e (t) — eV x u±(t)\ 

e^O Ae<±t<T ^ 

+ \\xu £ {t) - xu 6 ± (t)\\ L2 "j ^—-^ 0, 

where u± G C (R; S) are i/ie solutions to 

1 1 

ieStuI + -£ 2 Au| = F(x)n| ; u^, t=0 = ( 7 



L 2 



and ^j- are given by Proposition 



This result, as well as the analysis in can be viewed as a nonlinear analog to 
a result due to F. Nicr. In 42 , the author studies the problem 

ied t u e + \e 2 ^u e = V(x)u E + U Q u e ; uf t=0 = (^) , 

where [/ is a short range potential. The potential V in that case is bounded as 
well as all its derivatives. Under suitable assumptions, the influence of U occurs 
near t = and is localized near the origin, while only the value V(0) of V at the 
origin is relevant in this regime. For times e <C \t\ < T*, the situation is different: 
the potential U becomes negligible, while V dictates the propagation. Like in the 
nonlinear setting, the transition between these two regimes is measured by the 
scattering operator associated to U. 

The assumption a > tro(n) makes the nonlinear term short range. With our 
scaling for the nonlinearity, this perturbation is relevant only near the focus, where 
the potential is negligible, while the opposite occurs for e <C \t\ < T . 

The assumption on the potential V in Theorem 14.31 is such that this result 
is complementary to the analysis in j^J. It excludes the phenomenon of "total" 
refocusing. Indeed, the rays of geometric optics will all meet at one point for a time 
t =/= if and only if 5j = +1 for all j and the Wj's are pairwise rationally dependent. 
As discussed in [S] , this is the only case where the nonlinear term could be relevant 
at leading order, past the initial neighborhood of size e. 

Remark. An explicit change of variable makes it possible to state the theorem 
when </?(-) is replaced by ip( x ~ x ° )e Ix ^°/ £ (the phase factor e 1 *'^ / 6 is actually not 
relevant, since the only assumption we make is ip € S), see |12| . 

Remark. We did not mention the case a = 1 > na. Due to the lack of regularity of 
the nonlinearity F{z) — \z\ 2a z, this regime can be analyzed more easily when the 
power nonlinearity is replaced by a Hartree-type nonlinearity. In 1 1 1 j , we study the 
asymptotics behavior as e — > of the solution u e to 

ied t u" + ^e 2 Au^^u' + e a (\x\-^\u'\ 2 )u" ; uf t=0 = / , 

with 7 > 0, a > 1 and x £ M" (n > 2). In the same spirit as above, the expected 
critical values for the parameters are a > 1 or a = 1 to measure the relevance 
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of the nonlincarity outside of the focus, and a > 7 or a = 7 to decide whether 
the nonlincarity plays an important role at the caustic or not. When a = 7 > 1, 
results similar to Theorem 14.21 are proved. When a = 1 > 7, then as expected, 
the nonlincarity is relevant only away from the focus (and is measured by a slowly 
varying phase shift); the caustic crossing is described by the Maslov-Keller index. 
When a = 7 = 1, both nonlinear effects mentioned above are expected, and the 
underlying scattering operator would be a long range scattering operator. This case 
is treated only through a formal computation. 

Theorems 14 . 21 and 14.31 show that when the potential V is of the form 1*3.1*1 . the 

interactions between the linear dynamics and the nonlinear effects can be under- 
stood quite precisely. We mention two important limitations. First, we do not 
treat super-critical cases, such as A e = e a with a < na. The difficulty to treat such 
cases is not specific to the presence of the potential V, however. Second, a natural 
question is: what if V is a general smooth sub-quadratic potential? Some answers 
are given in |12| . As we already pointed out, the compatibility between linear and 
nonlinear analysis in [*3 El relies essentially in Lemma 14.11 In ["""""], we ask: for 
which potentials V can an operator of the form p(i)e*/ £ V x {e^/ £ ■) (compare with 
Lemma f4. II (ii)) commute with the linear part of i|1.2|l (point (i) in Lemma 1*4 . 1 f> ? 
Easy computations show that this is possible only if V is a second order polynomial 
(and when <f> solves the eikonal equation l"2.13|l ). For initial data of the form 14.2*1 . 
an interesting candidate as a substitute to J £ (which saves the day away from the 

focus in Theorems 14 . 21 and 14 . 3|) would be the Heisenberg observable U v {t)-U v {—t). 

1 2 

Note that it coincides with J e in the case V(x) = if 5 is taken as a new time 
origin. Obviously, it satisfies the commutation property (i) in Lemma 14. II It also 
yields a weighted Gagliardo-Nirenberg of the form stated in Lemma [4.11 (ii)', at 
least for \t\ < 8 small (the same as in Proposition 12.51) . In Lemma [4.11 the only 
interest of point (ii) is to imply (ii)' and (ii)", so only the action of this Heisenberg 
observable on nonlinearities of the form G(|z| 2 )z remains to be understood. This 
issue seems to be connected to Egorov theorem; we have no answer to provide. 



Guided by the semi-classical analysis described above, we can understand the 
role of some potentials on the finite time blow-up phenomenon. In the linear semi- 
classical analysis, it is well known that the energy is carried by the bicharacteristic 
curves. On the other hand, blow-up for (|2.1() is when the L 2 -norm of V x u be- 
comes infinite, while the L 2 -norm of u is constant; heuristically, the energy of u 
concentrates "somewhere" . One can expect that if in the linear case, bicharacter- 
istics meet, then the associated potential may encourage finite time blow-up in the 
nonlinear (focusing) case; the two effects cumulate. On the contrary, if the bichar- 
acteristics spread out, then the potential may compete with the attractivity of the 
nonlinearity. We give several illustrations that provide a rigorous support for these 
ideas. 

The first examples follow from the same heuristic, when one plays with the 
phase of the initial data instead of a potential in the equation. Consider the linear 
equation 



5. Changing finite time blow-up 




u £ (0,x) = f(x)e 




, b 61. 



The Hamilton flow is given by 



t = 1 



£ = + = 



X\ t =Q = X 



£|t=o = bx , 
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and the classical trajectories are: x(t) — xq{1 + bt). In particular x( = j^) = for all 
xq G R; if b < 0, then rays meet at the origin in the future, while if b > 0, they met 
at the origin in the past (see Figure EJ. Following the above heuristic discussion, 



t M t 




b>0 



Figure 1. Geometry of rays: case of quadratic oscillations. 

one may expect that quadratic initial oscillations alter the blow-up phenomenon in 
the nonlinear case. This is proved in |15| : let u solve (12.1(1 with uq e £, and suppose 
that blow-up may occur, that is A < and a > —. Denote u b the solution of ((2.1(1 

with initial datum Uq(x) = uo(x)e lb ^~ . Then for b ^> 1, u b is global in the future 
(t > 0). On the other hand, if Eq(uo) < 0, then we know from Proposition ^. 31 that 
u blows up at time T > 0, say. It is proved in |15| that if b < 0, then Ub blows up at 
time T b < = j^. For b < , the blow up phenomenon occurs sooner for u b than for 
it. In the conformal case a = —, the critical values for b can be explicitly related to 
the blow-up time T (see [T5l[Tlp . 
Consider now the equation 

ied t u £ + ^ 2 A« £ = (E ■ x)u e ; u E (0, x) = f(x)e ib1 ^ , b € R. 

The Hamilton flow is given by 

i = 1 ; x = C ; £ = -E ; f = ; x\ t=0 = x ; £|t=o = bx , 

and the classical trajectories are: x(t) = xq(1 + bt) — irE. Rays meet at x = —-^E 
at time This is the same phenomenon as above, shifted in space. Recall that 
we saw in Section 12.3.11 that the introduction of a linear potential E ■ x does not 
change the time of blow-up, but only shifts the set where this occurs. Here again, 
intuition and results meet. 

The last two cases we consider are isotropic harmonic potential and isotropic 
repulsive harmonic potential, with no initial rapid oscillation: 

1 \x\ 2 
ied t u £ ± + 2 £2Au ± = ±Uj2L ^~ u ± ; u ±(°> x ) = f( x ) ■ 

We have x(t) = Xoh(t), where h is given by 1(3.3(1 . Thus, x + (t) — xocos(ujt), 
and X-(t) = xo cosh(wt) . This is illustrated in Figure [21 Note the analogy with 
Figure ^ We already saw in Section 12.3.21 that the introduction of an isotropic 



NLS WITH POTENTIAL 



IS 




Figure 2. Geometry of rays: isotropic quadratic potential. 



harmonic potential may anticipate the blow-up time when a = — , just like quadratic 
oscillations in 15 . Similarly, the isotropic repulsive harmonic potential delays, or 
prevents, finite time blow-up. 

To complete the picture, we have to study the case A < 0, a > ~. 

Theorem 5.1 (0|B|). Let u G E, A < 0, a > \ and a < ^ if n > 3. Let u± 
be the solutions of i|1.2[l with e — 1 and V±(x) — ±w 2 -^-. 

1. If Ev + (uq) < ^uj 2 \\xuo\\ 2 L 2, then u + blows up at time 

2. If the initial datum uo satisfies 

,2 



< — 



5l|v«o|li> 



(T 



— ^Ilwolli^ < -yI|mo|I22, 



i/ien U- blows up in finite time, in the future or in the past. 
3. // the initial datum u$ satisfies 

,2 



|l|Vu ||| 2 



A 



G - 



1 



I w 



2cr+2 



i L 2„+2 < — HfmoIIl* 



Im / uqx ■ V x Uq 



then U- blows up in finite time, in the future and in the past. 
4. There exists mi > such that for any u> > u>i, the solution U- 



is global in time. 



Remark. Notice that in the first point, finite time blow-up occurs for a range of 
positive values of the energy (it was known that if Ey + < 0, then finite time blow- 
up occurs, see Jl]). This is in sharp contrast with the case of H2.1JI . where zero 
energy solutions may be global. When o~ = \ the solitary wave e lt R(x), where R 
is the ground state, solves l|2.1|l . is global in time and has zero energy. Note that 
the condition in the first point also reads Eq(uo) < (but Eq is not the energy 
associated to that equation!). 

The proof of points 2 and 3 relies on the Zakharov-Glassey method, just like 
Proposition ^. 31 so we shall not discuss it, and refer to The proofs of points 1 
and 4 rely on two conservation laws, which are more precise than the conservation 
of energy, and can be viewed as analogs to the pseudo-conformal conservation law 
()2.6(l . These laws have a geometric meaning, since they "follow" the propagation. 
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Lemma 5.2. Let uq G £, and a < — §^ i/n > 3. Lef it± &e the solutions of (|1.2I) 

I |2 

wii/i e = 1 and V± (x) — ±uj 2 . Introduce 

A 



£ 2 (f) 



2 

II 2 



El(t) 



(t)u+\\ 2 L 2 

(*)«H 
(*)«- 

( 

ff-(t)u_||£ 2 



cos 2 (^)||u+(i)| 



2o-+2 

£2<r+2 



2 

L 2 



■sin 2 (wi)||u+(t)||*£? a , 



■cosh 2 (w*)ll«-(*)ll£&, 
A^sinh 2 M)||,i_(t)|| 2 ^ 2 2 , 



where J+ stands for J when V(x) 
Ey ± . We have: 

dE]_ 



+uJ 2 ^-, and so on. Note that E\{t) + E^(t) 



(5.1) 
(5.2) 



UJ 



A 



dt 
dEl 
dt 



2ct + 

Uj\ 

2a + 



■{jig - 2)sin(2wi)||u+(t)|| 



£2<y + 2 



; (2 - ncr) smh(2u)t)\\u-(t)f L cr ++ 2 



Remark. These algebraic results can be proved in a classical way. Notice that one 
turn V+ into V- by replacing uj with iuj (and vice versa). 



The first point of Theorem l5 . ll follows easily. Assume that Ey + (uq) < 



\XUq 



2 

II 2 ' 



and suppose that u + exists in S up to time 
from (|5. ip . if cr > ^ and A < 0, 

dE, 1 



Then E\{Q) = E (u ) < 0, and 



dt 



< , Vi G 



IT ' 

°<2^ 



This implies 



E 



i fZ 
+ V2u 



2w/ 



< 0. 



But from the definition of E]_, 



Ei(^ 











XU+ 


\2uj) 









L 2 



so this leads to a contradiction (unless u + = 0, which means that uq = 0). There- 
fore, u + does not remains in S up to time It is easy to conclude that this is so 



because there exists T < such that 



lim 



\V x u+(t)\\ L 3 = oo . 



Notice that the larger uj, the sooner the blow-up. Note that the sufficient condition 
to have finite time blow-up does not depend on the value of uj > though. 

Even if the conditions 2 and 3 in Theorem 15. II become void as uj — » +oo, this 
does not mean that the last point is true. In the conformal case a — — , the 
last point is explicit, as we saw in 

Section EDO up to a characterization of finite 
time blow-up in this case where the potential is unbounded from below. Roughly 
speaking, the characterization of global existence is the same as in Theorem 12.11 
The energy Ey_ is the sum of three terms, and the term corresponding to the 
nonlinearity is controlled by the gradient and the L 2 norm of the solution, which 
is conserved. Thus, if the gradient remains bounded in L 2 , then each term of the 
energy remains bounded. Now the determinant of the matrix in (|3.5|l is constant, 
equal to — 1. Finally, to prove global existence, it is sufficient to check that the L 2 
norm of J_ (t)u remains bounded. 
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We briefly sketch the proof of the last point of Theorem 15. 11 and refer to [S] for 
details. The first step consists in noticing that the usual method to prove local 
existence still works for <|1.2|) with e = 1 and V(x) = — ^|a;| 2 , and does not "see" 
the parameter uj > 0. Indeed, as we already noticed, Mehler's formula yields the 
same dispersive estimate as in the case uj — 0, hence the same Strichartz estimates. 
Moreover, the operator J_ plays a role analog to that of in the case of (|2.1|1 : it 
commutes with the linear part of the equation, acts on the nonlinearity like a deriv- 
ative, and yields weighted Gagliardo-Nirenberg estimates with a weight uniformly 
bounded in w > (see Lemma [4.11 (ii)', and recall that coshx > 1). Therefore, 
there exists to > independent of uj > such that v.- € CQ — 2to, 2to[; S) and A(t)u 
is bounded in L°°([— to, to] : L 2 ) uniformly in w > for any A S {Id, J_, Now 
integrate (|5.2|) between to and t > to- Since A(2 — no) > 0, 

El(t) < Ei(t ) + Coj [ sinh(2ws)||w_(s)||^+ + 2 2 ds 
J to 

<C(C ) + Cu / 4 ^g^l|J-( S )-ll^- 
J t cosh{ujs) na 

The constants in the last estimate do not depend on to. Define 

y(t) := sup ||J_(a)«_||J a . 

t <s<t 

We have 

E 1 (t) < C(C ) + Cy{t)™' 2 f Sil ^ (2 ^ s) ds < C(C ) + = ? C : ry(i)" CT/2 . 

~ ^ ; yw 7 to cosh(ws) nCT ~ v ; cosh(^ ) n,T ~ 2 

We finally have 

y(t)<c(Co) + msh{ ^_M tr /2 . 

For cr > ^, we conclude by a bootstrap argument, since the constants and to do 
not depend on u > 0, and 

cosh(ujto) — > +oo . 

W — > + oo 

This yields a uniform bound for || J_(t)u-\\L 2 , and proves global existence. 

This proof relies on the evolution law H5.2J) , which seems to be bound to the case 
of isotropic potentials. A possible question is to ask whether a similar result holds 
when V is of the form i|3.1|) . with, say, 8% = — 1 (recall that the V is non-negative, 
things are rather well understood). An answer is given in |10| : 

Theorem 5.3 ([El). Take e = 1. Let n > 2, A e K, cr > f loiiA cr < 77=2 «/n > 3, 
anrf uo € S. Let V 6e 0/ £/ie /orm Ij3.1|l ™£/i <5i = —1, and denote 

u)± = max{ujj ; 5j = ±1} (lj + = if there is no 5j = +1) . 

Then there exists A = A(n, cr, |A|, ||uo||s) smc/i i/iai /or 

2cr 2 

(5.3) oj_ > A(l + w+) + — — (1 + oj+) ln(l + uj + ) , 

2 — (n — 2)er 

the solution u to (jl.2(l is global in time, u G C(M;S). 

The statement can be summarized as follows: if the repulsive force is sufficiently 
strong compared to other effects (linear confinement is overcome if oj_ 3> lo + , 
nonlinear effects are overcome if w_ ^> 1), then the solution is global. The strategy 
of the proof is as follows. First, in the same spirit as in [jj], we analyze the local 
existence result, to bound from below the local existence time, in term of the 
parameters Uj. Then, we notice that we obtain a time at which u is defined and 
small, and for which therefore the nonlinearity is not too strong. We consider the 
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solution of the linear equation f l|1.2() with A = 0) that coincides with u at that time. 
A continuity argument shows that u cannot move away too much from this linear 
solution. Since the linear solution is global, so is u. Note that the nonlinear term 
in uj + in (|5.3|) is zero when to + = 0, that is when there is no confinement. 

6. More on global existence 

In the last point of Theorem l5.ll and in Theorem 15. 31 we saw that if we consider 
a quadratic potential which is sufficiently repulsive, then the solution to l|1.2(l is 
global. In that case, we even have scattering: the solution to (|1.2|l is asymptotically 
linear, as time becomes infinite. 

It is proved in 8 that is the nonlinearity is defocusing (A > 0), and the potential 
V is the isotropic repulsive harmonic potential, then the solution u of 1|1.2|) is global, 
u S C(R; E). This follows from (JOJ- We then have: 

Proposition 6.1 ([HUHl)- Let e = 1, \,a> 0, with a<^ifn>3. 

1. Assume V{x) = -^-\x\ 2 . 

• For every it_ G £ ; there exists a unique uq G £ such that the maximal 
solution u G C(R;£) to l|1.2fl satisfies 

||CVH)u(t)-u_|| E — » o. 

t— * — oo 

• For every uq € £, £/iere exists a unique u + G £ smc/i ttcrf i/ie maximal 
solution u G C(M; E) to 1|1.2[) satisfies 

\\U v (-t)u(t)-u + \\ s — > 0. 

2. Suppose that V is of the form <Si = —1. 

• For every w_ G E, </iere ea;isi T finite and a unique u G C(] — oo,T]; E) ("1 
F 9 (] - oo,F];F 2ct + 2 (R")), where q = sokiion to (JT^ suc/i </iai 

||tV(-*)u(*)-u_|| s — » 0. 

• LetuQ G E. T/ien taking A larger in (|5.3[l if necessary, there exists a unique 
u + G E swc/i i/iai i/ie maximal solution u G C(R; E) to H1.2|l satisfies 

||tV(-*Mt)-u+|la^ 0. 

t— >+oo 

Notice that unlike in Proposition 12.41 there is no additional assumption on ct, 
simply a > 0: all power- like nonlinearities are short range for the Hamiltonian Hy 
under our hypotheses. 

The key ingredient of the proof consists in noticing that in the two cases, global 
existence follows from the boundedness of J(t)u in F°°(R;F 2 ). The weighted 
Gagliardo-Nirenberg inequality of Lemma f4. II yields exponential decay for u, and 
any positive power of an exponentially decreasing functions is integrable at infinity 

The restriction u G C(] — oo,F]; E) in the first part of point 2 comes from the 
fact that in general, we cannot prove that even if A > 0, the solution u is defined 
globally in time. 
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